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ABSTRACT 

Context. Solar cycles vary in their amplitude and shape. There are several empirical relations between various parameters 
that link cycle's shape and amplitude, the foremost the Waldmeier relations. 

Aims. The solar cycle is believed to be a result of the solar dynamo action, therefore these relations require an explanation 
in the framework of this theory, which we aim to present here. 

Methods. We related the cycle-to-cycle variability of solar activity to fluctuations of solar dynamo drivers and primarily 
to fluctuations in the parameter responsible for the recovery of the poloidal magnetic field from the toroidal one. To 
be specific, we developed a model in the framework of the mean-field dynamo based on the differential rotation and 
a-effect. 

Results. We demonstrate that the mean-field dynamo model, which is based on a realistic rotation profile and on 
nonlinearity that is associated with the magnetic helicity balance, reproduces both qualitatively and quantitatively the 
Waldmeier relations observed in sunspot data since 1750. The model also reproduces more or less successfully other 
relations between the parameters under discussion, in particular, the link between odd and even cycles (Gnevyshev-Ohl 
rule). 

Conclusions. We conclude that the contemporary solar dynamo theory provides a way to explain the cycle-to-cycle 
variability of solar activity as recorded in sunspots. We discuss the importance of the model for stellar activity cycles 
which, as known from the data of the Mount Wilson HK project, which measures the Ca H and K line index for other 
stars, demonstrate the cycle-to-cycle variability similar to solar cycles. 
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1. Introduction 

Solar activity has a periodic nature, but the cycle ampli- 
tude and shape vary from one cycle to the other. This chal- 
lenges the prognostic abilities of solar activity models. The 
sunspot activity can be quantified by using various trac- 
ers derived from observations. These tracers show an inter- 
relation among each other. The indices characterizing the 
tracers can be employed t o pred ict the future evolution of 
solar activity. Waldmeier ( 1935 ) first suggested this option 



(an inverse correlation between the length of the ascending 
phase of a cycle and the peak sun spot number of that cycle) 
and applied it fWaldmeier"1936) to predict the subsequent 
cycle. Later, other relations of this kind wer e proposed and 



called Wa ldmeier relations (see fo r a review Vitinsky et al 



m 



(see discussion, 
|etal.,2008i Karak & C: 



e^ 



1986 and Hathaway et al. 2002). The nature of physical 
processes underlying t 



le Waldmeier relations is not clear 



Cameron & Schiissler 2008; Dikpati 

^ !nQudh uri 2011). We note, that these 

statistical properties of the magnetic activity also exist for 
other tracers related to the sunspot activity (e.g., sunspot 
group number and ar ea, see [Vitinsky et a l. 1986; Hathawayj 



et al. 

other 



| 2002| [Karak fc Choudhuri| ] 

ind of solar and stellar activity indices, e.g., for the 



2011[ ), and even for the 
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Ca II index ( Soon et al.|1994 ). The Waldmeier relations are 
conside red as a valuable test for dy namo models (Karak fc[ 
|Choudhuri|[2011|[Pipin fc Kosovichev[[2011a[ ). 

Clarifying the physics underlying the Waldmeier rela- 
tions is particularly attractive to support the prognostic 
abilities concerning solar cycle. It is believed that the cyclic 
solar activity is driven by a dynamo, i.e. a mechanism that 
transforms the kinetic energy of hydrodynamic motions into 
a magnetic one. Many modern solar dynamo models (see, 
e.g., ,Stix.2002} assume that the toroidal magnetic field that 
emerges on the surface and forms sunspots is generated near 
the bottom of the convection zone, in the tachocline or just 
beneath it in a convection overshoot layer. This kind of dy- 
namo can be approximated by the Parker dynamo waves 
(Parker 1993|). The direction of the dynamo wave propa- 
gation in the framework o f the aCl-dynamo is defined by 
the Parker-Yoshimura rule ('Parker'1955; 'Yoshimur a|1975 ), 
according to which the wave propagates along iso- surfaces 
of the angular velocity. The propagation can be affected by 
the turbulent transport (associated with the mean drift of 
magnetic activity in the turbulent media by means of the 
turbulent mecha nisms), by the anisotropic turbulent dif- 
f usivity (Kitchatinov[[2002 ), and by meridional circulation 
( [Yoshimura ' 1975[ Choudhuri et al.] [l995). An alternative 
to the Parker's surface dynamo waves is the distributed 
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dynamo with subsurface shear (e.g., Brandenburg 2QQ5| ), 
where the dynamo wave propagates along the radius in th e 
main part of the solar convection zone ( Kitchatinov||2QQ2| ) . 
Near-surface activity is determined by the subsurface shear. 
Another popular option is the flux-transport dyn amo (e.g., 
Choudhuri et al.||1995| [Dikpati fc Charbonneau| [l 999 ) . 

In the context of dynamo theory, the Waldmeier rela- 
tions can be explained by invoking physical mechanisms of 
the solar magnetic fleld generation and a mechanism that 
drives variations of the amplitude and shape of t he activity 
cycle. For example, Pipin & Kosovichev (2011a), hereafter 
PKll, showed that variations of the a-eflect amplitude may 
explain the correlation between the cycle rise rate and the 
cycle amplitude and other t ypes of the Wald meier relations 
as well. It was sues; ested (Choudhuri 1992; Hoyng| p93| ) 
that the fluctuations of the a-effect (associated with ki- 
netic helicity fluctuations) are likely to be one of the natural 
sources of the cycle parameter variations. 

In addition to the statistical relations between the cycle 
parameters within a separate cycle there are correlations 
relating the parameters in subsequent cycles, for example, 
the odd-even cycle and the last cycle period-amplitude ef- 
fects. These effects are closely related to the memory of 
the dynamo processes and to the strength of the saturation 
processes, which damp deviations of the cycle parameters 
from the refere nce state characterizing the cycle attractor 
( Ossendrijver & Hoyng 1996; Ossendrijver et al. 1996). 



It was argued (Choudhuri 1992 



|Hoyng|[l993l ), that a 

dozen percent is a reasonable estimate for tne noise com- 
ponent of the a-effect. Previous calculations (see the above 
cited papers) showed that a straightforward application of 
the idea with the vortex turnover time and the vortex size 
as the correlation time and length for the a-fluctuations 
needs fluctuations much stronger than the mean a. On the 
other hand, the results of direc t numerical simulations (e.g, 
Brandenburg fc Sokoloff||2QQ2 ) and results of current helic- 
ity (related to a) o bservations in solar active regions (e.g., 
Zhang et al. 2010) suggest that the correlation time for 
^-fluctuations can be comparable to the cycle length and 
the correlation length comparable to the extent of the lat- 
itudinal belts. Using t hese results. Moss et al. (2008) and 
Usoskin et al. (2009b) showed that an a- noise on the or- 
der of few dozen percents is sufficient to explain the Grand 
minima of solar activity. The aim of this paper is to exam- 
ine the result of a-fluctuations in the statistical properties 
of the solar cycle including the Waldmeier relations and the 
odd-even cycle effect. 

We chose a particular model for the solar cycle in which 
a-fluctuations are introduced. Of course, it is impractical 
to try all available models to learn which one is better to 
obtain the relations under discussion, but we select below 
the model among a relative wide choice of the models that 
gave bet ter results in the preliminary simulations (jPipin^^ 
Sokoloff|[2QTT] ). 



2. Basic equations 

2.1. Dynamo model 

The dynamo mod el employed in this paper ha s been de- 
scribed in detail in Pipin & Kosovichev ( 2011a|b ). We study 



the standard mean-fleld induction equation in perfectly 
conductive media: 



dB 

~dt 



V X (£:+U X B) , 



where £^ = u x b is the mean electromotive force, with u, b 
being the turbulent fluctuating velocity and magnetic fleld, 
respectively; U is the mean velocity (differential rotation), 
and the axisymmetric magnetic fleld is 



B esB + V X 



r smt 



where is the polar angle. The expres sion for the mean 
electromotive force vector S is given by Pipin (2008). It is 
expressed as follows: 



£i = {(^ij + (1 + ^7) 7ij) 5 - (1 + ^77) r]ijkVjBk. 



(1) 



Tensor a^j represents the alpha effect. It includes the hy- 
drodynamic and magnetic helicity contributions. 



Co.(l^^a)^a{P) Sin' ^C^lf 



(M) 
• cy ■ ■ 



(2) 



^if^X 1 where the 



where the hydrodynamical part of the a-effect is deflned 
by OLlf\ ^Q!,r7,7 deflnes the noise, ^jj^ {P) is the quenching 

1^1 

function, where /3 = \: _ , u is the convective RMS 



velocity. The reader can find expressions for the quenching 

function, ifjoc and a-^^ in appendix. 

Contribution of the small-scale magnetic helicity x = 

a-b (a is a fluctuating vector-potential of the magnetic fleld) 

to the Qf-effect is deflned as a)^^ 

''J 

coefflcient c\^^ depends on the turbulent properties of the 
medium and on the parameter characterizing the influence 
of the Coriolis force on convection. Expression for c\^^ is 
the same as in PKll and is given in appendix, as well. Other 
parts of Eq.([T]) represent the effects of turbulent pumping, 
7ij, and turbulent diffusion, rfijk. We give their expressions 
in appendix. 

The nonlinear feedback of the large-scale magnetic fleld 
to the a-effect is described as a combination of an "alge- 
braic" quenching by the function ifjcy {j3)^ and a dynamical 
quenching due to the magnetic helicity conservation con- 
straint. The magnetic helicity, x ? subject to a conserva- 
tion law, is desc ribed by the following equation (Klee orin] 
fc Rogachevskii|fl999^ iSubramanian fc Brandenburg,,2004p : 



dx 
dt 



-2 (5.B) 



R^Tc 



V • (% VX) : 



(3) 



where Tc is a typical convection turnover time. Parameter 
controls the helicity dissipation rate without specifying 
the nature of the loss. It seems reasonable that the helic- 
ity dissipation is most efflcient in the near surface layers 
because of a strong decrease of Tc toward the surface. The 
last term in Eq. ([s]) describes the diffusive flux of the mag- 
netic hehcity (Mitra et al.J|2010|). 

We u sed t he solar convection zone model computed 
by Stix ( |2002[ ), in which the mixing- length is deflned as 

I = OiMLT |A*^^^| , where A*^^^ = Vlogp is the pressure 
variation scale, and aMLT = 2. The turbulent diffusivity is 



parametrized in the form, tjt = Cr^r]^\ where 77^^ = 
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Table 1. Parameters of the dynamo models 



Model 




Ca/Rx 


Bo 


noise 


a 


2Dla 
2Dlry 
2DlaL 


10-^ 

-/- 
-/- 


0.03/200 
-/- 
-/- 


800 
-/- 
-/- 


exp(^c.) 


0.2 
0.04 

logcr(^c.) 



is the characteristic mixing-length turbulent diffusivity, i is 
the typical correlation length of turbulent flows, and is 
a constant to control the efficiency of the large-scale mag- 
netic field dragged by the turbulent flows. Currently, this 
parameter cannot be introduced into the mean-field theory 
in a consistent way. 

In this paper we use = 0.05. The fairly low tur- 
bulent diffusivity both due to low <C 1 and due to 
quenching of the turbulent difffusivity coefficient for the 
fast-rotation regime in the depth of the solar convection 
zone, where = 2Qorc ^ 1, provides the correct value of 
the cycle period in the model. Note that in the fast-rotation 
regime the turbulent magnetic diffusivity is dominated by 
the anisotropic component of the diffusivity tensor along 
the rotation axis. This component is growing in the inter- 
mediate range variations of 1^*. Currently, the problem with 

<C 1 has no satisfactory resolution within the dynamo 
theory. In the model we used an analytical fit to of the dif- 
ferential rotation profile proposed bylAntia et al. (1998). 
It was given in our earlier papers (see Fi g, la in ^Pipin fc, 
Kosovichev||2Qlla| [Pipin fc Sok5bff||2QTT| ). 



We matched the potential field outside and the per- 
fect conductivity at the bottom boundary with the stan- 
dard boundary conditions. For magnetic helicity, similar to 
Guerrero et al.| (|2Q1Q) we put V-rX = at the bottom of 
the domain and = at the top of the convection zone. 

The parameters of the model are summarized in the 
Table [ij where r]^/r]T is the ratio between the turbulent 
magnetic helicity diffusivity and the turbulent magnetic dif- 
fusivities; the parameter controls the helicity dissipation 
rate; Bq is a typical strength of the toroidal magnetic field 
controlling the sunspots number parameter in the 2D mod- 
els; the column "noise" defines the fluctuating parameter 
and a is the standard deviation of the Gaussian noise in 
the model. The lognormal noise in the model 2DlaL was 
symbolically denoted as exp(^Q,). 

The left panel in Figure [l] shows a typical time-latitude 
diagram in the model 2D la for the toroidal magnetic field 
evolution averaged over the subsurface layers 0.9 — 0.99Rq 
and the radial magnetic at the top of the integration do- 
main. The right panel shows the time-radius diagram for 
the toroidal an poloidal components of the large-scale mag- 
netic field evolution at 30° latitude. Note that the geometry 
of the poloidal magnetic field is drawn by iso-contours of 
the toroidal vector-potential. 



u 




Fig. 2. Typical realization of the fluctuating part of the a- 
effect (left panel) and its probability distribution function 
(solid line, right panel). There we show the PDF for the 
lognormal fluctuations as well (dashed line). 



cause a'^^^ ~ u' and tjt ^ u''^)- To examine the long- 
term dynamics of the model with regard to the specific 
statistical distribution of the noise we included the results 
for a mod el wit h the lognormal distribution of (see. 
Subsection 3.2.3). In this case the parameters of the log- 
normal distribution were computed from the corresponding 
Gaussian distribution. Random numbers were generated 
using the Numerical Recipes Fortran subroutine ''gasdev". 
The amplitudes of the fluctuations were restricted to 2(j. 
A realization of the lognormal fluctuations was prepared 
before the run. We took the input parameters from the 
realization of the Gaussian distribution fluctuations, which 
were computed with taking into account the 2a cut-off. The 
fluctuation renewal time was constant and equal to the pe- 
riod of the cycle in the model. Figure |2] shows a typical 
realization for the Gaussian fluctuations of the a-effect, its 
probability distribution function (PDF), and the PDF of 
the generated lognormal fluctuations. 

We considered the amplitude of the standard deviation 
of the a-coefficient as an input parameter. We chose this 
quite arbitrary based on the following crude estimation. 
The total number of cyclones in the solar co nvection zone 

2 X 10^, see, e.e;., Miesch et al. 



may be of the order of = :^ x iU^, see, e.g. ^ 

(2008]), who found that the solar convection vertical vortic- 
ity spatial spectrum is flat and has a maximum at about 
£ ^ 140, where £ is the spatial wavenumber. If a ^ O.lu^^ 
then the fraction of the magnitude fluctuations in 

velocity corresponds to about 20% of the fluctuations in 
a. We stress that this very crude estimate needs a support 
from detailed numerical modeling and observations, which 
are, however, obviously out of the scope of this paper. The 
Gaussian fluctuations were cut off on the 2cr-level to exclude 
a possible effect of rare and very strong fluctuations, which 
are hardly associated with Waldmeier relations. It can be 
important in principle for Grand minima statistics. We did 
not address the latter possibility systematically, but tested 
the lognormal fluctuations to study the impact of the PDF 
tail on the long-term solar cycle variations. 



2.2. Noise model 

In Eq.([T]), the noise, ^^,77,7, contributes to the hydrody- 
namic part of the a-eSect (see, Eq.([2|), to the turbulent 
diffusion, and to the turbulent pumping. Following Usoskin] 
et al. (2009b) the model employs the long-term Gaussian 
fluctuating ^a,r] of the low amplitude with RMS deviation 
given in the Table 1 (last column). It is expected from 
general consideration that, for the low amplitude fluctu- 
ations, is an order of magnitude smaller than ^q, (be- 



2.3. The sunspot cycle model and the Waldmeier relations 

Here, we define the Waldmeier relations as a set of the mean 
properties of the sunspot cycle including relations between 
the rise rate of the cycle and its amplitude and the relation 
"period-amplitude" . In the original form the Waldmeier re- 
lation reads as a link between the period of a cycle and 
amplitude of the subsequent cycle. Other relations like this 
(rise time and amplitude; rise time and decay time) are 
sometimes referred to Waldmeier relations, as well. These 
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3050 3055 3060 3065 3070 3075 3080 3050 3055 3060 3065 3070 3075 3080 

TIME, [YR] TIME, [YR] 



Fig. 1. Typical time-latitude and time-radius (at the 30° latitude) diagrams of the toroidal field (gray scale), the radial 
field (contours at left panel) and the poloidal magnetic field (which is drawn by the contours of the vector-potential at 
the right panel) evolution in the 2D la model (see Table 1). The toroidal field averaged over the subsurface layers in the 
range of 0.9 — O.99i?0, the radial field is taken at the top of the convection zone. 



rela tions were considered i n PKll and in our previous pa- 
per ( Pipin fc Sokoloff|2011 ). The ratio of the decay and rise 
rates is known as the shape of the cycle. Amother relation 
present ed in literature is the so -called Gnevyshev-Ohl rule 
(e.g. Charbonneau et al. 2007'), which provides a positive 
correlation between the amplitudes or intensities of 2A^th 
and 2N + 1th cycles. 

The amplitude of a cycle is defined as the difference 
between the maximum sunspot number and the sunspot 
number in the preceding minimum. The latter can differ 
from zero because of the overlap of subsequent cycles. The 
cycle period is equal to the time between the subsequent 
minima. The rise time of a cycle is defined by the differ- 
ence between the moment of the cycle maximum and the 
moment of the preceding minimum. The rise rate is defined 
as the ratio between the difference of the sunspot numbers 
at maximum and minimum of the cycle and the rise time 
of the cycle. There is similar definition for the decay rate 
of the cycle. 

Following to PKll we relate the sunspot number with 
the toroidal magnetic fields in the near-surface layer. In 
PKll the instant sunspot number was defined using the 
maximum of the toroidal magnetic field strength, which 
was taken over all latitudes and averaged over the surface 
layer from 0.9 to O.99i?0. Then, this value was related to the 
sunspot number (SSN) via the three-halves law proposed by 
Bracewell (1988) to relate SSN in the cycle with a "linear" 
sinusoidal part of the SSN variation. The relation between 
the toroidal field and SSN, which was introduced in PKll 
has an undesirable property of giving non-zero SSN at the 
minimum of a cycle and strong variations of the minima 
amplitudes with variations of the a-effect parameters (see 
Fig. 6 there). 

Here, we exa mined another possibility (also see [Pipin 
fc SokolQff||2011 ). We assumed that sunspots are produced 
from the toroidal magnetic fields by means of the nonlinear 
instability, and avoided to consider the instability in details. 
To model the sunspot number W produced by the dynamo 
we used the following ansatz: 



W (t) Cw {Bma^)sLexp 



Bo 



(4) 



where (^max)^/, is the maximum of the toroidal magnetic 
field strength over latitudes averaged over the subsurface 
layers in the range of 0.9 — 0.99Rq; Bq is a typical strength 
of the toroidal magnetic field sufficient to produce the 



sunspot, hereafter we put Bq = 800G; Cw is the param- 
eter to calibrate the modeled sunspot number re lative to 
observat i ons. H ereafter we put Cw = 1- Results by [Pipin fc 



Sokoloff[ ( [2QTT] ) showed that similar to PKll the Waldmeier 
relations can be reproduced with the Wolf number defini- 
tion (Eq|4|. The exponential dependence in Eq.Q yields 
the Waldmeier relations at smaller variations of the a-effect 
compared to those in PKll, where the Bracewell law was 
used. Also, we find that the mean-field dynamo model with 
relation Q reproduces the long-term variation of the cycle 
much better than in the case of the three-halves law. 



2.4. Observational data set 

Although the series of group sunspot numb ers covers 400 



years smce 1610 AD ([Hoyt fcSchatten|1998'), givmg a mea- 
sure of the temporal variability of solar activity, parame- 
ters of the solar cycle such as its total length and ascend- 
ing/descending phases are not reliably known for the earlier 
times. Solar cycle parameters can be more or less reliably 
evaluated since 1750 or, with some caveats, after the end of 
the Maunder minimum in 1712 (Usoskin 2008). However, 
even in this case an uncertainty related to the po tentially 
lost solar cycle in t he last decade of 18-century (e.g. [Usoskin 



et al.[[200 3 

observed 



2009a) still exists. The Sun was amazingly well 
uring the Ma under minimum , especially in 



its 



second half (Ribes fc Nesme-Ribes [199 3| , but the solar cy- 
cle was suppressed below the threshold for sunspot forma- 
tion, which led to unclear dynamo manifestation s. Cycles 
before the Maunder minimum are not well known ( [Vaquero 
et al.'2011) and their shapes cannot be obtained. Therefore, 
only the period of 1750-2009 AD, which includes 23 full so- 
lar cycles in the official numeration, can be analyzed here. 
Statistical properties of the long-term variations of the so- 
lar cycle can be esti mated on the base of the reconstructed 
data set proposed by Usoskin et al. (2004) and Solanki et al. 
( 2004) . 

There are several synthetic series that present solar 
cyclic variability for the times before the beginning of the 
sunspot series. They are based on a fit of a prescribed math- 
ematical model to fragmentary non-systematic qualitative 
prox y data of naked-eye sunsp ot or auroral observation 



(e.g., 

ries 



Schove 1955 



LSpC 

Nagovitsyn] 1997 ) . These synthetic se- 
do not pretend to be quantitative reconstructions of 



solar activity and cannot be used to analyze of solar cycle 
parameters, which are explicitly prescribed in the model 
rather than reconstructed. 
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Although sunspot activity is greatly suppressed during 
Grand minima, the solar dynamo continues to operate at a 
redu ced level. For example, an analysis of sunspot and au- 



rora ( Kfivsky fc Pejml|l988 ) data during the Maunder min- 
imuni suggests that the dominan t periodicity was shifted 



Table 2. The parameters of the cycle given by the dynamo 
medel 2D la and by the monthly smoothed actual sunspot 
number (SSN). 



from t he 11 years to 20-22 years ( Silverman]! 199 2 Usoskin 
et al.p QQl). Data of the cosmogenic isotope ^"^C also con- 
firm longer cycles du ring the Maunder minimum ( Peristykh 
fc Damoii[l 998; Mi yahara et al.|2QQ6b ). A similar lengthen- 
ing of the solar cycle during a Grand minimum has been ob- 
served, using the ^^C data, also for the Sporer minimum a t 
the turn of the 15-16-th centuries (Miya hara et al.| |2QQ6a). 
However, the parameters of individual solar cycles cannot 
be recovered for the Grand minima periods, only the sta- 
tistical features. 

Taking into account all the above information, we com- 
pared simulations with the monthly smoothed sunspot 
number (SSN) data set from the SIDC (Solar Influence 
Data Center), which starts at 1750. The data set was 
additionally smoothed by means of the Wiener filter. To 
compute the wavelet spectra of the sun spot dat a se t, we 
used the data set provided by |Hoyt fc Schatten| ([1998| and 



Solanki et aL] ( |2QQ4| 



3. Results 

We performed long-term simulations for the time interval 
of about 10^ years (i.e. the time-sp an of the longest re con- 
struction of the solar cycle history ( Solanki et al.||2QQ4 ) us- 
ing our basic model. To compare the results obtained with 
other dynamo models see Pipin fc Sokoloff (2011). The sim- 
ulated time data set for W{t) is shown in Fig. (top panel). 
It shows events with the extended period of the high- and 
low magnetic activity. 



3.1. The Waldmeier relations and the odd-even cycle effect 

First of all, we divided the set into separate cycles to 
look for the Waldmeier relations (rise rate to amplitude 
and period to amplitude) and the odd-even cycle effect 
(Gnevyshev-Ohl rule). The results are shown in the bot- 
tom panel of Fig. |3j which depicts spreads of simulated 
points over corresponding planes, and their linear fits (solid 
lines) . The dashed lines show the correlations computed on 
the base of the actual sunspot data. The results for the rise 
rate to decay rate relation is similar to that for the rise rate 
to amp litude relation and can be found in [Pipin fc Sokoloff] 
(2011). Data concerning the linear fits shown in Fig. are 
given in Table [2j where the first four rows contain infor- 
mation for the mean and variance (standard deviation) for 
the parameters of the sunspot cycles in the different data 
sets. The shape of the cycle is defined as the ratio between 
the decay rate and the rise rate of the cycle. The last five 
rows show linear fits with the mean-square error bar and 
the correlation coefficient for the Waldmeier relations and 
for the odd-even cycle effect, (I) marks the effect that is 
calculated on the base of the SN integrated over the cycle 
and (II) marks the effect that is calculated on the base of 
the cycle amplitudes.. 

We see from the Fig. |3] that the model reproduces the 
Waldmeier relations and the Gnevyshev-Ohl rule reason- 
ably. Note that the dispersion of both the simulated and 
observational data from the linear fit of the rise rate to 
amplitude (as well as that one for rise rate to decay rate 





2Dla 


SSN 


Period 


11.07±1.08 


11.01±1.12 


Amplitude 


103.3±40.5 


108.2±38.1 


Rise Time 


4.06±.77 


4.32±1.07 


Shape 


.59±0.15 


.69±0.31 


Rise Rate - 


3.3x+14.8±6, 


2.9x±33.2±8.9, 


Amplitude 


0.98=1= 0.001 


0.97± 0.01 


Period - 


-17.4X+277.6 


-23.6x±368.5 


Amplitude 


±27.5, 


±28.0, 




-0.54=1= 0.02 


-0.68± 0.12 


Rise Time - 


-43.1X+259.4 


-26.7x±224. 


Amplitude 


±24.2, 


±25., 




-0.67± 0.02 


-0.76± 0.1 


Odd - Even(I) 


0.68x±155.±136., 


0.35x±235±145, 




0.67± 0.03 


0.33 ± 0.3 


Odd - Even(II) 


0.58x±35.6±26.6, 


0.35x±62.5±32.3, 




0.58± 0.03 


0.42± 0.28 



Pipin fc Sokoloff (2011)) are much lower than that for the 
the period to amplitude relation and Gnevyshev-Ohl rule. 
We composed (see Fig. 2, bottom right) a relation of the 
rise time to amplitude (i.e. using the quantity inverse to 
the rise rate) to learn that the dispersion looks more or less 
like that for the relation perio d-amplitude. As discussed by 
Cameron fc Schiissler (2008), a relation rise time to am- 
plitude has higher dispersion than that for the rise rate 
to amplitude. We conclude that the quality of fitting sub- 
stantially depends on the presentation method chosen to 
illustrate a relation. 

The results concerning the even-odd effect are shown in 
Fig[4j Note that the original formulation of the Gnevyshev- 
Ohl rule is calculated on the basis of the sunspot number 
(SN) integrated over the cycle (Fig. [4] right), while the left 
panel of the figure shows the effect calculated for the cycle 
amplitude. We see from the figure that the model repro- 
duces observations in both cases more or less reasonable, 
however, the slope substantially depends on the definition 
of the Gnevyshev-Ohl rule. 

It is expected that the strength of the sunspot cycle 
depends on the strength of the poloidal field of t he Sun in 
the precedin g solar minimum. Following this idea, |Schatten 
et al. (1978) suggested using the strength of the Sun's po- 
lar field for the cycle prediction. Recently, this idea was 
exploite d in t he Babcock-Leighton type model, see the re- 
view by Hathaway (2009 ) . 

Figure |5j illustrates the phase relation between the am- 
plitude of the sunspot cycle and the strength of the dipole 
component of the dynamo-generated magnetic field. There 
we show the backward and forward correlation between 
these parameters of the model, taking the strength of the 
dipole component at the cycle minimum. The strength of 
the dipole component refers to the surface, and it was cal- 
culated as the first coefficient in the spectral decomposition 
of the magnetic potential A. The backward correlation has 
the correlation coefficient 0.86 ± 0.01 and approximation 
145 Ax — 74. ± 13. The forward correlation between the cy- 
cle amplitude and the strength of the dipole component 
of the dynamo-generated magnetic field at the subsequent 
minimum has the correlation coefficient 0.64 ±0.01 and ap- 
proximation 0.004x ± 0.8 ± 0.12. This relation has a higher 
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Fig. 3. Top, the simulated W{t) in our dynamo model. Bottom, the Waldmeier relations for the model; squares show 
data for individual cycles, while the solid line gives the correlation, the dashed line shows these relations obtained from 
the actual SSN data. 
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Fig. 4. Odd-even cycle effect (Gnevyshev-Ohl rule): the left panel shows the effect that is calibrated on the base of the 
cycle amplitude, the right panel shows the effect that is calibrated on the basis of the SN integrated over the cycle. 
Squares show data for individual cycles, while the solid line gives the correlation, the dashed line shows these relations 
obtained from the actual SSN data. 



dispersion than the results the fluctuation of the a-effect. 
For comparison we added several points obtaine d from the 



WSO polar magnetic field observations ( [Svalgaard et al. 



1978[ |Hoeksemanl995^ 



3.2. Other perspectives 

The main aim of this paper is to demonstrate that fluctua- 
tions of the a-coefficient provide an option to explain short- 
term dynamics of solar activity cycle such as Waldmeier 
and similar relations. We note, on one hand, that this idea 
can be useful to explain more long-term dynamics and, on 
the other hand, that a is far from being a unique transport 
coefficient in dynamo equations, which can be noisy. These 



6 



Pipin et al.: Variations of solar cycle profile in a solar dynamo 




AXIAL DIPOLE, [G] AXIAL DIPOLE, [G] AMPLITUDE, [SSN] 

Fig. 5. Left, the phase relation between the cycle amplitude and the strength of the dipole component of the dynamo 
generated magnetic field at the surface. Middle, correlation between dipole component at the cycle minimum and the 
amplitude of the subsequent cycle. Right, correlation between the amplitude of the cycle and the dipole component at 
the subsequent cycle minimum. The circles mark the results of the WSO observations and the SSN data. 



noisy transport coefficients can obviously contribute the ac- 
tivity cycle dynamics. Of course, a detailed investigation of 
these options is far beyond the scope of this paper, but we 
present below some exploratory results in these directions. 



3.2.1. 7^-fluctuations vs ^-fluctuations 

Obviously, fluctuations of the a-coefficient affect the solar 
activity evolution together with fluctuations of other dy- 
namo governing parameters. Comparing the relative role 
of various fluctuations in the solar cycle variations is be- 
yond the scope of this paper and we restrict presentation 
by comparison of the effect of r]- and a-fluctuations only. 

Choudhuril (|1992|) addressed this point and suggested 



the following relation for the fractional growth rate F of per- 
turbations in the dynamo as a result of the perturbation of 
the Qf-effect and the turbulent diffusion (in our notations): 



P 



p2 



where P is the period of the cycle and Td is the typical dif- 
fusive time of the dynamo. This relation gives a hint that 
the fluctuations of the turbulent diffusivity are relatively 
more significant for the dynamo perturbation provided the 
values S,a/<^ and £,etci/T] are comparable. The point is, how- 
ever, that a ^ u while r] v? ^ where u is the turbulent 
velocity. We performed simulations and present their re- 
sults in Fig. [g] for our model with a-fluctuations of 20% 
(Fig. |6j left) and 4% /^-fluctuations (Fig. |6] right), which 
corresponds to the comparable fluctuations in u. We see in 
Fig. |6] that the a-fluctuation looks more pronounced than 
/^-fluctuations. 

We examined the influence of the pumping effect fluctu- 
ations on the cycle variations. The parameters of the model 
are the same as for 2D la, except that we let the a-effect 
be constant and has the same characteristics as ^q, in 
2D la. It is found that for a model with fluctuations the 
cycle variations look very similar to the model 2D1?7. The 
model demonstrates that the cycle amplitude variations are 
approximately half of those for the 2D la model. Variations 
of the period are also quite small, like in the model 2D1?7. 




Fig. 6. Cycles distribution for the ISOOyr data set from 
2Dla model (^c^-noise), left and 2Dl7^ model (^^-noise), 
right. 



3.2.2. Resonance effects 

Our base model exploits the memory time of a-fluctuations 
equal to the nominal cycle length. A natural worry is that 
a resonance could participate in the Waldmeier relations 
simulated while the correlation time of the a-coefficient in 
solar convection can be different from the cycle length, thus 
avoiding resonances. Note that the resonance effects for dy- 
namo waves is almost not addressed in scientific literature 



(Gilman fc Dikpati 2011). We varied the correlation time 
and calculated the cycle amplitude variance (Fig.[7|). Some 
peaks are visible in this figure, which may indicate reso- 
nance effects, however, the renovation times to which they 
are attached vary from one run to the next. On the other 
hand, the results given in Figure [sj right) suggest that the 
resonance effects may depend on the phase synchronization 
between the fluctuations of the a-effect and the cycle vari- 
ations. Therefore, we may anticipate that the fluctuations 
on the descending phase of the solar cycle are more effective 
than those on the rise phase of the sunspot cycle. Bearing 
in mind the distributed character of the dynamo model, we 
conclude that the resonance phenomena that po ssibly play 
a role here ne ed to be addressed separately (cf. Gilman & 
Dikpati||2QllD . 



3.2.3. Long-term dynamics 

We move from the dynamics activity cycles, viz. the 
timescale of several dozens cycles where the Waldmeier re- 
lations and the even-odd relation is applicable, to the much 
longer term history of the solar cycle on timescales of up to 
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Fig. 7. Dependence of the cycle amplitude on the noise re- 
newal time. Three runs with various renovation times are 
shown with solid, long-dashed, and short-dashed lines re- 
spectively. 



10^ years. Here we cannot discuss such fine details as the 
Waldmeier relations because the available data does not 
trace the cycle shape. Th e available isotopic d ata (we refer 
to the reconstruction by Solan ki et al.| l2QQ4) only traces 
the evolution of the cycle- averaged quantities. Because the 
present dynamo model was adopted (via tuning the param- 
eters Ca^Rx arid T]^) to reproduce short-term dynamics, 
we expect that the long-term dynamics will be reproduced 
not as well. 

Figure [S] shows the global wavelet power spectra for the 
sunspot data and for the model data sets. The reader can 
find the definition of the concept of wavelet power spectra 
and its discussion in the context of the solar activity studies 
in, e.g., Frick et al. (1997a). The data were processed with 
the standard MATLAB wavelet toolbox using the standard 
Morlet wavelet basis. The global wavelet power spectra 
were obtained by integrating the spectra in the time do- 
main (see Eq.(lO) in the cited paper). Each spectrum was 
normalized relative to its maximum magnitude. To illus- 
trate the role of the a-effect fiuctuation statistics we show 
the results for the model with the log- normal noise ^q,. The 
mean and variance of the log-normal noise corresponds 
to the mean and variance of the Gaussian noise in the 
model 2Dla. The short-term scales spectra look qualita- 
tively similar in all three sets. The principal difference is the 
ratio between the spectrum amplitude for the basic cycle 
(at 11 years) and the amplitude of the second maximum at 
the period ^ 200 years. This ratio is greater in the sunspot 
data set. All three data sets show the long-term variations 
with periods of about 100 years, which is usually identified 
with the Gleisberg cycle. 

The dynamics on the scale of millennia looks similar for 
all three data sets, see Figure [s] (right). The model with 
the log-normal noise does not show the ordered long-term 
variations on this time scale, while the model 2D la and the 
reconstruction data show evidence for the variations with 
period about 6000 years. It is unclear, however, if this result 
is statistically stable. 

An important statistical property of the dynamo is the 
occurrence pro bability of hi gh- and low-activity episodes. 
Following Solanki et al. ( 2004J , we defined the high-activity 
episode as having the average SN > 50 (the minimal aver- 
aged SN was higher than 50, min((5'A^)) > 50) and similar 
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Fig. 9. Probability of the high- and low- activity episodes 
for the given duration (see definition in the text). The thin 
lines show the results for the high-activity episodes, and 
the bold lines the low-activity episodes. The solid lines (thin 
and bold) show the results for 2Dla model, the dotted lines 
show the results for the 2DlaL model and the dashed lines 
show the same for t he reconstruction data set provided by 
Solanki et all (|2004l). 



for the low-activity episode occurs when max((5'A/')) < 50. 
Then, we counted the number of episodes for each time 
scale with the high and low activity episodes and computed 
the probability distributions as a function of the time scale. 
The results are shown in Figure [9) We find that the dynamo 
models show a somewhat higher probability for the high- 
activity episodes than the reconstruction data set. Their 
probability profiles looks similar in all three cases and show 
a significant drop in the pass from decadal to the centen- 
nial time scale. The probability of a high-activity episode 
to occur decreases exponentially with time. The probabil- 
ity of the opposite event, i.e., a minimum with the average 
SSN below 50, increases accordingly. Note that the end of 
a high-activity episode does not necessarily imply a low- 
activity epoch. It could be a local minimum with a duration 
of less than 20 years (the 11-year interval was used for av- 
eraging). A similar behavior is found for the probability of 
the low-activity episode in the dynamo models that show 
a significant drop of the probability profiles around half- 
millennium. The reconstruction data set is very different in 
this aspect. It seems that it is not possible to explain all 
the basic properties of the sunspot cycle variations as the 
result of fiuctuations of the a-effect. 

Can one predict that the high-activity episode will be 
ended with the low-activity episode of the given length? 
Figure [To] shows the results for the low-activity episodes of 
30 and 50 years. The given estimates can be biased because 
of the unpresentable statistics for these events. 



4. Discussion and Conclusions 

We have studied the impact of low-amplitude Gaussian 
fiuctuations of the a-effect on the statistical properties of 
the magnetic dynamo cycle, such as the Waldmeier rela- 
tions and the Gnevyshev-Ohl rules. The dynamo model in- 
cludes long-term fiuctuations of the a-effect and employs 
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Fig. 8. Wavelet spectra of the simulated and observational sunspot number data sets. The left panel corresponds to short 
time scales of up to 300 yr. The sol id line shows the r esults for the 2Dla model, the dashed line is computed on the 
base of the data set provided by the Hoyt & Schatten (1998) reconstruction and the dash-dotted line shows the results 
for the 2DlaL model based on the log- normal fluctuations of the a-effect. The right panel is simil ar to the left one, bu t 
for longer time scales. The dashed line is computed on the basis of the reconstruction provided by Solanki et al. (2004). 
The spectra were normalized relative to their maxima for clarity. 




60 

TIME, [YR] 

Fig. 10. Probability that the high- activity episode will be 
ended by low-activity episode of a given length, 30 years 
length - bold lines and 50 years length - thin lines. The 
solid lines show the results for 2D la model, the dotted 
lines show the results for the 2DlaL model based on the 
log-normal fluctuations of the a-effect, and the dashed lines 
show t he res ults for the reconstruction provided by |SolankT| 
eFaIl(l2004l). 



two types of a nonlinear feedback of the mean-field on 
the a-effect, including algebraic quenching and dynamic 
quenching due to the magnetic helicity generation. The gen- 
eral properties of the dynamo, such as the direction of the 
toroidal magnetic field drift, the polar magnetic field sign 
reversal at the maximum of a cycle, etc., are consistent with 
observations. 

Our model does not include the meridional circulation 
effect, which is advocated by the Babcock- Leighton and the 
fiux-transport type dynamo models (e.g., Choudhuri et al. 



1995| [Dikpati fc Charbonneau | |1999 
It was shown that a part of the Wa. 



Dikpati et al.||2008 ). 

dmeier relations can 



be possibly explained by a specially tuned fiux-transport 
model that considers fluctuations of th e meridional circula- 
tion speed (Karak & Choudhuri 2011). Still, observational 
constraints on the distribution of the meridional circulation 
inside the convection zone are not very strong because we 
have measurements for the surface. The angular momen- 
tum balance in mean-field models supports for the circula- 
tion pattern, which has a deep circulation stagnation point 
and a strong concentration of the velocity speed towards 
the bottom and the top boundaries of the s olar convection 



zone (e.g., Kitchatinov & Olemskoy 



|201lD. Yet, most of 
IKarak 



& ChoudhurjJ 



the flux-transport models (including | 
2011) use a very different circulation pattern. Following to 
this argumentation, we postpone a more complete study of 
the effects of the meridional circulation fluctuations to the 
future. 



We showed, con firmi ng the previous findings of Pipin 



KosovichevI ([2011a|) and [Pipin fc Sokoloff| ( |201lD , that vari- 
ations of the a-effect amplitude result in variations of the 
cycle amplitude and period. Taking into account random 
fluctuations of the a-effect, we calculated statistical prop- 
erties relating the cycle amplitude, the cycle shape, the rise 
time, etc., on the basis of the simulated SN data set cov- 
ering period of more than 10000 years. Our results agree 
well with observations for the Waldmeier relations and the 
Gnevyshev-Ohl rules. 

From the qualitative point of view these results were an- 
ticipated from the earlier ana lysis of the helicity fluc tuation 
effect in the dynamo give n by Choudhuri' ('1992') and'Hoyng^ 
1993) (see, also Ossendrijver fc Hoyng 1996; Ossendrijver 



^ir[|l996| [Moss et al.||20Q8| [Usoskin et al.||2Q09bp . Our 

results presented in Figure [5[ about the correlation of the 
polar dipole field and the cycle amplitude and as the results 
for the Gnevyshev-Ohl rules suggest that the Waldmeier re- 
lations can be understood by considering the general prop- 
erties of the magnetic field generation processes, which are 
involved in the dynamo. 
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Our model shows a good correlation (with low variance) 
between the strength of the polar dipole magnetic field in 
the cycle minimum and the amplitude of the subsequent 
cycle. This results from the deterministic process of the 
toroidal magnetic field generation by the differential rota- 
tion from the large-scale poloidal magnetic field . This cor- 
relat ion is often used for the cycle prediction (Hathaway 



300 



|2QQ9| ) by Babcock-Leighton type dynamo models and it is 
for tne first time demonstrated in the mean- field dynamo. 
The rise rate of the sunspot cycle depends on the differential 
rotation and the amplitude of the poloidal field. Therefore, 
the correlation between the rise rate and amplitude of the 
cycle is a derivative property and is a consequence of the 
link between the polar dipole magnetic field in the cycle 
minimum and the strength of the toroidal field in the sub- 
sequent cycle. 

Fu rther more, follow i ng the general idea of [Zaslavsky] 
( |1978| ) (cf, |Hoyng||1993[ |Charbonneau et alT]|2QQ7| ), we can 
interpret the Gnevyshev-Ohl rule as evidence that the so- 
lar cycle is a nonlinear self-excited oscillation that tends to 
preserve the property of the attractor under random pertur- 
bations. The amplitude and phase of the subsequent cycles 
are related by the so-called Zaslavsky map. The strength of 
the link between the parameters of the subsequent cycles 
is controlled by the fluctuation amplitude and by the per- 
turbation's decrement. The latter strongly depends on the 
nonlinear mechanisms involved in the dynamo. To exam- 
ine this idea we made additional simulations with a lower 
helicity dissipation rate (high R^) and found that the cor- 
relation coefficients between the parameters of the subse- 
quent cycle increase with increasing R^. Therefore the link 
between the odd and even cycles, and the period to ampli- 
tude correlation in subsequent cycles can be considered as 
evidence for the fluctuation impact on the dynamo and ev- 
idence for nonlinear damping of these perturbations in the 
dynamo. This conclusion needs to be investigated in more 
detail especially by comparing the results of the a-effect 
and meridional circulation fluctuations. 

Long-term variations of the magnetic cycle in the dy- 
namo can be induced in different ways. Two main mecha- 
nisms can be identified: nonlinear deterministic chaos and 
an effect of fluctuations of the turbulent parameters in- 
volved in the dynamo process. Generally, we anticipate that 
statistical properties of long-term cycle variations are de- 
pended on the force that drives the long-term variations. 
We examined weakly nonlinear models with the amplitude 
of the a-effect close to the threshold. In our models, the typ- 
ical ratio between the energy of the large-scale toroidal field 
and the kinetic energy of convective flows does not exceed 
0.3. As a result, the chaotic regime in the model is not as 
evident as the impact of the a-effect fluctuations. Figure |9] 
illustrates the difficulty to obtain extended episodes of low 
magnetic activity in this case, w hile these episodes are com- 
mon in the solar dynamo (Solanki et al. |2004J . To amplify 
the chaotic regime in the model, we have tried additional 
possibilities and included an angular momentum balance 
into the dynamo problem to take into account the nonlin- 
ear feedback of the magnetic field on th e differential rota - 
tion. The model was described earlier by<Pipin ( |1999 2004). 
Figure 11 shows the simulated SN for the model involving 
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Fig. 11. Time series of the simulated sunspot number for 
the extended 2Dla model involving the fluctuations of the 
alpha effect and the magnetic feedback on the differential 
rotation. 



occurrence of lowactivity episodes as that in the reconstruc- 
tion data set. 

It is natural to expect that at least stellar magnetic cy- 
cles of solar-like stars should demonstrate a variability sim- 
ilar to the solar one, including relations comparable with 
the Waldmeier relations. Available observations of stellar 
activity provide some hints that support this expectation. 
Stellar activity data of the Mount Wilson HK project mea- 
suring the Ca H and K line index for other stars (Baliunas| 



the nonlinear effect of the magnetic field on the angular mo- 
mentum balance in the solar convection zone. This model 
shows higher intermittency in the cycle variations than that 
in Figure [3) and indeed it has a similar probability of the 



et al.|1995 ) are available for two activity cycles. The wavelet 
analysis of the data for several stars (Frick et al. 1997) 
demonstrated that the subsequent cycles for a given star 
can differ in their cycle amplitudes. We note that a moni- 
toring of stellar activity of solar-like stars to obtain relations 
similar to the Waldmeier ones could establish our prognos- 
tic abilities of solar activity based on these relations much 
better. 

Summarizing the results of the paper, we conclude that 
the mean-field solar dynamo theory provides a way to 
explain the cycle-to-cycle variability of solar activity as 
recorded in sunspots. The results given in the literature 
and the results obtained in the paper suggest that the 
Waldmeier relations can be explained invoking very differ- 
ent kinds of dynamo models. More work is necessary to 
study the relations between the statistical properties of the 
dynamo cycle, and the dynamo mechanisms involved in the 
magnetic activity will help to obtain more insight into the 
processes operating in the stellar and solar dynamo. 
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Appendix A: Appendix 

Here, we describe the contributions of the mean- 
electromotive force th at are i nvolv ed in Eq.([l]). The basic 
formulation is given in |Pipin| ( [20081 ) (PQ8). In th IS paper we 



(H) 

reformulate tensor ol\ j \ which represents the hydrodynam- 
ical part of the a-effect, by using Eq.(23) from FOB in the 
following form, 



o^' = % {^VT {Ao' (e • A(^)) + /I") (e • A(«))) } (A.l) 



The contribution of magnetic helicity x = ^'b (a is a fluc- 
tuating vector magnetic field potential) to the a-effect is 



defined as a, 



(M) 



C'jf^X^ where 



(x) 



(A.2) 



The turbulent pumping, 7^ j, is also part of the mean elec- 
tromotive force in Eq.(23)(P08). Here we rewrite it in a 
more traditional form (cf, e.g., ), 



7.i = 3r?r {/i"^A(f) + (e ■ A('')) e„} e, 



(A.3) 



(a) 



The effect of turbulent diffusivity, which is anisotropic be- 
cause of the Coriolis force, is given by 



Vijk 



= 3vt { (2/1"^ - fi''^) e,,k - 2/l"^e,e„£„,fc} . (A.4) 



Functions f^i'flij depend on the Coriolis number ft* = 
2rcfto and the typical convective turnover time in the 
mixing- length approximation, Tc = i/u' . They can be found 
in F08. The turbulent diffusivity is parametrized in the 



form, r]T 



Cr^r]^\ where r]^^ 



is the characteristic 



mixing-length turbulent diffusivity, vi! is the RMS convec- 
tive velocity, i is the mixing length, is a constant to 
control the i ntensity of tu rbulent mixing. The other quan- 
tities in Eqs.( A.l|A.3|A.4 ) are A*^^^ = Vlogp is the density 
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stratification scale, A*^^^ = Vlog (j]^"^^ is the scale of tur- 
bulent diffusivity, e = fi/ \Q\ is a unit vector along the 
axis of rotation. Equations ( A.1|A.3|A.4 ) take into account 
the influence of the fluctuating small-scale magnetic fields, 
which can be present in the background turbulence and 
stem from the small-scale dynamo. In our paper, the pa- 

rameter e = — which measures the ratio between the 

magnetic and kinetic energies of fluctuations in the back- 
ground turbulence, is assumed to be equal to 1. This corre- 
sponds to the energy equipartition. The quenching function 
of the hydrodynamical part of a-effect is defined by 

Note in the notation of P08 tpa = ~3/4(/)g^\ and /3 = 
1^1 



Below we give the functions of the Coriolis number 
defining the dependence of the turbulent transport genera- 
tion and diffusivities on the angular velocity: 

/!"' = (3 + + 10.- 5) 

- ((8e + 5)^^*2 + 30e- 15)) , 

- ((4£ + 5)Si-= + 15{e-l))), 
r = -:^('3 ((3.-11, fi-+5e- 21) 



- (4 (£-3)1]*^ + 156 -63)), 

- {5{e + l)n*^ + 12e-Q)) , 



- {{2s + l)n*^ +3S-3)) , 
/^i^ = -^f3(l.- + l)(l.- + 2.-l--*-^ 
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